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{Sj ' We perform a lattice study, in the quenched approximation, of dynamical mass generation in 

a system of relativistic (Dirac) fermions, coupled to an Abelian gauge field in (2 + l)-dimensions, 
in the presence of an external (constant) magnetic field, perpendicular to the spatial planes. It 
is shown that a strong magnetic field catalyzes chiral symmetry breaking, in agreement with 
results in the continuum. The role of the higher-Landau poles in inducing a critical temperature 
above which the phenomenon disappears is pointed out. We also discuss the implications of this 
model on the opening of a gap in doped antiferromagnetic superconductors. 
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1 Introduction 



In (2 + l)-dimensional gauge theories with relativistic fermions chiral symmetry can be 
defined only in the case of an even number of massless fermion species, 2N, in the so- 
called 4x4 reducible Dirac formalism M. Fermion mass terms break the chiral symmetry 
U(2N) — > U(N) x U(N). To preserve parity one needs to have masses with opposite sign 
between the fermion species. This is the only type of mass terms that is allowed to be 
generated dynamically in vector-like theories, as a result of energetics arguments . 

The problem of chiral symmetry breaking in (2 + l)-dimensional continuum QED 3 in 
the presence of external magnetic fields has been studied extensively in the recent liter- 
ature 0, £|, following analyses in four dimensional gauge theories |)J [J[ ||. According to 
these findings, a strong enough magnetic field acts as a strong catalyst of chiral symme- 
try breaking, which occurs even in the limit of weak gauge couplings. We remind the 
reader that, in the absence of such strong fields, dynamical chiral symmetry breaking (i.e. 
fermion mass generation) occurs only in supercritical theories i.e. for dimensionless 
gauge couplings that are stronger than some critical value g c Q. The effect of the strong 
magnetic field is, according to the analysis of |6| , to reduce g c to zero. This phenomenon 
was termed in || magnetic catalysis. 

This phenomenon may be important in the physics of high-temperature supercon- 
ductivity, as pointed out in ref. ||. Indeed, it is believed that high-temperature super- 
conductors may have some strong connection with relativistic (2 + l)-dimensional gauge 
theories 0, 10, [11], || • This belief is enforced recently, as a result of experiments indicating 



that the high-temperature superconductivity is actually due to the opening of a <i-wave 



gap in the hole (fermion) spectrum [12]. As is well known, c?-wave superconductivity is 
characterized by lines of nodes of the gap. In the gauge theory approach to doped an- 
tiferromagnets, such hole (charged) excitations about such nodes are described by Dirac 
fermions || p| , coupled to gauge fields that describe the magnetic interactions among 
the charged degrees of freedom, in a spin-charge separating formalism [ 13] . 

The models of ref. ||, || describe the dynamical opening of a fermion gap at those 
nodes. Moreover, the theoretical analysis of || indicated that the dynamical fermion gap, 
due to the statistical gauge interactions, is enhanced in the presence of strong external 
magnetic fields, as expected from the general theory of || |3J, |J . The analysis of [§] also 
indicated the existence of a critical temperature above which chiral symmetry is restored 
in the model. This phase transition occurs well within the superconducting phase of the 
model and pertains only to the fermion gap at the nodes of the ci-wave gap Q. For the 

1 In (2 + l)-dimensions the effective dimensionless gauge coupling is actually the inverse of the fermion 
flavour number, l/N. Thus supercriticality in the coupling means actually the existence of a criti- 
cal number of flavours iV c , below which dynamical chiral symmetry breaking is possible. Continuum 
Schwinger-Dyson analyses 0, in the large N approximation, have indicated that N c ~ 32 /n 2 . This 
result has also been argued to be correct on the basis of lattice analyses . 

2 We note that the problem of looking at charged femions in the presence of an external field is com- 
patible with the superconducting nature of the high-temperature superconductors. These materials are 
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model of ||, in the absence of external magnetic fields, the fermion gap at the nodes 
disappears at temperatures 0.1 K at most; in contrast, under the influence of a strong 
magnetic field 1 Tesla < B < 20 Tesla, the fermion gap at the nodes remains up to 
temperatures of order 30 K. 

It is worth pointing out that recent experiments have indicated the presence of such 
gaps at the nodes of a conventional c?-wave gap [13 ] . Such experiments involved magnetic 
fields of the above order of magnitude, and they showed clearly the existence of a gap up 
to temperatures that are in close agreement with the above theoretical results. Although 
there may be alternative explanations for the phenomenon, the above qualitative agree- 
ment between theory and experiment prompts further studies of the magnetic catalysis 
phenomenon. 

In addition to the above-mentioned motivation from Condensed Matter Physics, the 
study of the magnetic catalysis phenomenon in three-dimensional gauge theories is also 
useful from the point of view of high-temperature four-dimensional gauge theories in the 
Early Universe, or the physics of high-temperature quark-gluon plasma, in the presence of 
external fields, since both cases may be effectively described by three-dimensional gauge 
models. 

Whatever the motivation might be, it becomes clear that the problem of the formation 
of a fermion condensate in the presence of strong external fields is quite relevant for 
physical applications. So far, as far as we are aware of, the analysis was restricted to the 
continuum formalism ||, It is the purpose of this work to present a first, preliminary, 
analysis of the phenomenon on the lattice, in the quenched approximation in the case of 
(2 + l)-dimensional QED. 

We shall discuss here the formation of parity-invariant fermion condensates. In our 
analysis we shall trigger the phenomenon of chiral symmetry breaking by adding bare 
masses which preserve parity, and then extrapolate the results to the zero bare mass case. 
In this way we find only magnetically-induced condensates which preserve parity. Notice 
that in the presence of an external magnetic field, the theory ceases to be vector-like, and 
hence the theorem of no longer applies, and one could in principle have parity-violating 
fermion condensates. This phenomenon was conjectured by Laughlin [15| as happening 
in the case of the high-temperature superconductors, offering a possible explanation of 
the results of the experiments of [fI4}| . Our results in this work, though, indicate that 
parity- violating induced condensates might not be necessary for a possible explanation of 



the situation in the experiments of [14 



known to be strongly type II, which implies that the Meissner effect penetration depth in the supercon- 
ducting phase is quite large, a few thousands of Angstroms. 
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2 Results in the continuum from a novel perspective 

In this section we describe briefly the theoretical formalism in the continuum |3], |j. One 
follows essentially the method of Schwinger flEf , by looking at the coincidence limit of 
the fermion propagator (in configuration space), Lim x _> y S(x,y)\s, in the presence of a 
constant external field, B. We start first from the free-fermion case, i.e. the case where 
the fermions interact only with the external constant field. The external gauge potential 
is then given by: 

A„ = -Bx 2 5^ (1) 

and the Lagrangian is: 

L = -^(zY^/i - ieAf) - m)tf (2) 

where m is a parity conserving bare fermion mass, and the 7-matrices belong to the 
reducible 4x4 representation, appropriate for an even number of fermion species formal- 
ism \U 
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The induced fermion condensate is given by the coincidence limit of the fermion prop- 
agator 

< 0|1>$|0 > \ B = +Lim x ^ y trS(x,y)\ B 

where S(x,y) =< 0\T^{x)^{y)\0 > \ B (3) 



Following the proper time formalism of Schwinger jTJj], the propagator S(x,y)\s in 
the presence of a constant external magnetic field can be calculated exactly ||: 

S{x,y)\ B = exp (lej^Afdz^ S(x - y)\ B 

~ S W\b = I 8(7r t )3 /2 ex P[- ? (j + sm ')] x exph^x.^%] 

x[(m + y^C^x" - ^F™^ x (esBcot(eBs) - y^T**?)] (4) 

where is the Maxwell tensor corresponding to the external background gauge po- 
tential (0), and = rf v + ([F ext ] 2 y u [l - eBscot(eBs)}/ B 2 , and the line integral is 
calculated along a straight line. A useful expression, to be used in the following, is the 
Fourier transform Se(A:)|b in Euclidean space of S(x — y)\s 0: 



Se(^)\b — — i dsexp[—s(m 2 + k 2 + k\ 
Jo 



tanh(eBs) . 



eBs 1 

x{[— k^ + m — z(& 2 7i — ki^2)tanh{eBs)} x [1 — i^^tanh^Bs)}} (5) 

We define 

A < > \ B =< 0|tftf|0 > \ B - < 0|^|0 > \ B =o (6) 

3 Wc define for convenience the condensate with a plus sign, rather than the usual minus sign. 
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A straightforward calculation, then, yields in (2+l)-dimensional space-time the following 
result for A < ^ > \ B 0: 

A<m>\ B = J d 3 ktr[S E (k)\ B - S E (k)\ B=0 ] 

[/ " / d 3 k f ds[exp{-s(m 2 + kl + k\ tanh ^ Bs ) ^ _ exp {_ s ( m 2 + kl + k 2 )}} 
J Jo eBs 
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(27T) 

777,o R /*oo p — sm I 

/ ds^—rco^eSs) - (7) 
Jo V s e -o s 



2ttI 



In the formulae above we have used the notations: k\ = k\ + k\. The limit of the 
final expression as m — > is easily found to be We note here that the infinities in 
< > are already present for zero magnetic field, so by performing the subtraction to 
form A < > 1 73 we get an automatically finite result and there is no need to introduce 
explicit cut-offs in the integrals. 

We show, in figure 1, the outcome of the numerical calculation of (0) for three different 
masses. The magnetic field in the following figures will be represented by the ratio b = — . 
Q There is a striking result, depicted in this figure, namely that A < > \ B increases 
for decreasing mass. This may be understood through the remark that the response of 
the system to an external field is bigger for small masses. An important point is that the 
value of this quantity has a finite limit as m — > 0, the value referred to above. 

Further insight in the model, as well as the tools to investigate the finite tempera- 
ture issue, is offered by the expansion of the propagators in terms of the Landau level 
contributions. The relevant expansion reads: 



g ( 1)n D n (m,B,k) 



S E (k)\ B = -ie-rt V(-l)" " ' ' , (8) 



n=0 



k 2 + M 2 (B) 



with M 2 (B) = m 2 + 2eBn and TrD n (m, B, k) = 4m[L„(§) - L n _ x (2*L)]. 

We need the integrals over the momenta fci, k 2 , which are calculated using the formula: 

k 2 , Ik 2 2k 2 
dk x dk 2 e-^[L n {-^) - L„_i(-i)] = (-l) n 2neB. (9) 
eJD et) 

What is left is the sum over the k momentum. We consider the case of zero temper- 
ature (the case we have just treated), which means one should perform the integral over 
k . The result reads: 



eB emB ~ 1 eB mV2eB ^1 



m 



2 



< W > \ B = — + ^ — — = -— + —z Cb, ttb) ( 10 ) 

^ n 71 n~=i M n [B) 2tt 2tt 2 2eB 
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The denominator tt has been included here to facilitate the comparison with the lattice results: there 
the quantity b is defined as the ratio e g\ B ^ — , where eB\ max is the maximum magnetic field (in lattice 

units) that can be achieved on the lattice; this maximal value is tt. 
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Figure 1: < > versus — at zero temperature 



where ((z, q) is one of the Riemann ( functions, which is understood to be defined through 



appropriate analytic continuation |17|. The first term in the middle part of the equality 
fllOD is mass independent and coincides with the previous result in the massless limit, while 
the remaining part of the expression vanishes in this limit. Thus, for the zero-temperature 
case, we find that in three dimensional gauge theories a strong magnetic field (such that 
\eB\ is much larger than any other mass scale in the problem), may induce a fermion 
gap Q. 

It is worthy of examining some physically interesting limiting cases of ([TOD, keeping 
however eB finite. First, consider the case of weak magnetic fields (compared to the (bare) 
fermion mass m): eB << m 2 . I that case one may expand M" 1 ~ mT 1 (l — + . . .V 
The resulting sums can be regularized by means of the Riemann's second ( function, 

-§> C( 



CO*) = E 



oo 
n=l 



n 



using 0: C(0) 



1 



-1/12 (both being understood by 



5 Notice an important difference from the corresponding four-dimensional problem, where the corre- 
sponding fermion condensate reads: 



rrf) A 2 

< >~ ~\eB\ — (In— + <D(m )) 
Air'' m z 











(11) 



and tends to zero as the bare mass m — > , for a fixed (ultraviolet ) scale A. However, the quantum 
dynamics of the electromagnetic field may drastically change the situation if treated non-perturbatively [[| 
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means of appropriate analytic continuation). The result is: 



A<VV> U« m2 ~ + ... (12) 

12m z 7c 

with the . . . indicating subleading terms of higher order in eB/m 2 . Thus, for weak mag- 
netic fields the induced gap scales quadratically with the external field intensity, for mas- 
sive fermions. 

In a similar manner, for strong magnetic fields eB » m 2 , one finds 

A < « > \ eB>>m , = LJ. [l + ct-)^ - ppfl^) + ■ ■ J (13) 

where £(-) — —1-46, £(|) — +2.61.. Notice that despite the fact that the various terms in 
(P~3|) come with alternating signs, however there is no critical value of the field, consistent 
with the approximation m/y/eB « 1, at which the condensate in (^) changes sign. 

We now come to the model at finite temperature. The integral over fc should be 
replaced by a sum over ko = 2n(n + |)T. The T = result is modified to: 

- . eB , ,m, emB ™ 1 7/ M n (B), 
< „ > | T = -««*(-) + — E — (anM^), (14) 

thereby implying the absence of the condensate for any finite temperature if the bare 
infrared cut-off mass m, for finite constant B. 

In a similar fashion, as for the zero-temperature case, one may perform the summation 
over the Landau poles for certain limiting cases. The physically interesting case is the 
one for which one assumes a strong field eB » m 2 and a temperature regime T » m, 
but such that eB » T 2 . In such a case tanh(M n (B)/2T) ~ 1 for all n, whilst M n (B) ~ 
\/2eBn, to leading order. Hence, 

— eB m 1 .46mV eB 
A<«>| T ^_- ^— (15) 

using C(l/2) = —1.46. Notice the existence of a critical temperature, for non-zero m, 

T c ~ 1 (16) 

above which the condensate vanishes 0. The order of magnitude of the temperature is 
consistent with the approximations made in the derivation of (|T6"D, which suggests an 



6 To be precise, at T = T c the condensate changes sign. However, given that in three-dimensional 
gauge theories < WW >oc E 2 , with E the (magnetically-enhanced) fermion energy gap, such a change in 
sign is interpreted as implying that above T c thermal fluctuations dominate, and, hence, prevent chiral 
symmetry breaking. It is in this sense that we used above the terminology 'critical temperature'. 
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important role for the higher Landau levels at finite temperatures, [] namely that they 
prevent dynamical mass generation above a certain temperature. This is consistent with 
the fact that for B — > oo, where only the lowest Landau poles are retained, there is 
dynamical mass generation for all T < oo. In the context of a possible application of 
this phenomenon to high-temperature superconductors ||, we note that V eB scaling of 
a critical temperature with the external field intensity are reported in the experiments of 



14| . Notice however, that the existence of extra interactions, that may lead to dynamical 
mass generation for fermions, as in the SU(2) x Us(l) model of [IT], the scaling of the 
critical temperature with B may be different ||. 

We now proceed with a different regularization, based on the remark that all infinities 
in the model lie in the T = sector, so that the difference 

Tt . , ^ T . > eB 1 2emB ~ 1 1 

Trl + er tt n=1 M n {B) i + 



of the condensate at finite temperature, as given by (0), and its zero temperature limit 
should be finite. On the other hand, the zero temperature part has already been regular- 
ized above by the subtraction of the B = part and the finite result is given in equation 
(0). Thus, we may calculate A < > | T substituting for A < tyty > \ T=0 the regular- 
ized expression of (|7|) and for the finite difference the series of (|ITD- The result of these 
calculations is depicted in figure 2. 

We have used two masses, m = 0.20 and m = 0.05. We remark that the regularization 
scheme we have adopted seems to work quite well. We observe that (say for the larger 
mass, where the effects are more transparent) one may distinguish three temperature 
regimes, delimited by T\ = Ta ~ 0.04 and t<i = Ta ~ 0.17. Below T\ the condensate is the 
same as for zero temperature, thus if one uses lattices of time extent equal roughly to 
which equals about 25 in the m = 0.20 case, one should get zero temperature results. After 
that the condensate starts decreasing and vanishes at ti. Thus, for lattices with temporal 
extent of about ^ (about 6 in the m = 0.20 case) one should expect the vanishing of the 
condensate. For the smaller value of the mass, 0.05, the zero temperature value is larger, 
in agreement with the results in figure 1, but it decreases more steeply. The corresponding 
values for T\ and r<i are 0.01 and 0.145 respectively, yielding the corresponding temporal 
extents 100 and 7. Comparing the results for the two masses we observe that for small 
temperatures (corresponding to large enough lattices) the condensate for the smallest 
masses are the largest, as explained in figure 1; however, if the lattice is relatively small, 
corresponding to a high temperature, it is expected that the condensates for the "big" 
masses will take over. This could be used as a check that the lattice size used is large 
enough. 



7 Especially as far as the superconductivity phenomenon is concerned, there is a natural upper bound 
for the magnetic field, B c , above which superconductivity is destroyed. For realistic cases, the critical 
field in the high-temperature superconducting oxides is of order of 20 Tesla. For such fields the inclusion 
of the higher Landau poles seems necessary. 
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Figure 2: < > versus temperature at fixed magnetic field 



We have plotted the value of the critical temperature T c (denoted by r 2 above) versus 
the magnetic field in figure 3. 

We have used data for two masses and performed a square root fit: T c = K(m)Vb = 

0.52, so that 



K,(m) 
«(0.20) 



(eB). It turns out that «(0.05) = 0.46 and k(0.20) 



0.26, 



^ — 0.29. We observe that the fit is better for the small mass, which is quite reason- 
able, since the square root behaviour has been found in the limit of big magnetic field as 
compared to the mass squared. The big values of the magnetic field have been fitted, so 
the discrepancies occuring in for the bigger mass take place at small magnetic fields. The 
proportionality constants 0.26 and 0.29 given by this approach is equal to a very good 
approximation with the value \ yielded by the ( function regularization fllEp, applicable 
when the magnetic field is strong. We note that strong field means value of the ratio ^| 
much greater than 1. In the case we are considering, the value b = 0.1 corresponds to 
^§ ~ 8, for m=0.20 and to ^4 — 128, for m=0.05. Thus the two different approaches 
agree quite well and reproduce the magnetic field dependence found experimentally. 
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Figure 3: Critical temperature versus magnetic field for two masses: the upper data and 
fit correspond to m = 0.20, the lower data and fit to m = 0.05. 



3 The Lattice Model 

We now describe briefly the field-theoretical model associated to the problem of super- 
conductivity as formulated above. The relevant three-dimensional Lagrangian is: 

L = -\(fp 2 +^D^ - mm, (18) 

with the covariant derivative. The field strength corresponds to the "statistical" gauge 
field a^. The mass term has been introduced by hand to facilitate the lattice treatment 
and the limit m — > will be considered in the end. It has already been explained that the 
Clifford algebra is reducible in three dimensions and the model may also be represented 
by 2 x 2 matrices acting on two species of fermions with opposite mass parameters. The 
generators 1, 73, 75 and A = 27375 generate a global SU{2) x U{1) symmetry of the 
Lagrangian (if m is set to zero). We observe that one may write down two kinds of mass 
term: the parity conserving mass term mc$$ and the parity violating one: myf Af . It is 
easy to see that the parity conserving mass term breaks the global SU (2) symmetry, while 
the parity violating mass is an SU(2) singlet. The U(l) field acting on the fermions may 
consist of two parts, the electromagnetic 6^(1), generated by an external electromagnetic 
potential A^, and the aforementioned statistical Us(l), associated with a'^, responsible 
for the fractional statistics of the holon excitations of the doped antiferro magnet. 
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The model with an external homogeneous magnetic field has been examined in || in 
the two-component formalism for the fermions; it was found that the (parity violating) 
magnetic field gives rise to parity violating mass terms with opposite signs for the two 
fermion species. Thus, in the four-component formalism, the dynamically generated parity 
violating mass terms are expected to sum up to zero. On the other hand, the effect of 
mass generation will manifest itself in an enhancement of the parity conserving mass, as 
compared to the free fermion field case. 

We would like to study the above effects by lattice methods; this will free us from the 
need of the approximations used in || and will permit us to include the statistical gauge 
field, on top of the external electromagnetic field. The lattice action, corresponding to 
the continuum Lagrangian, may be written as: 

S = fa £(1 - TrUp) + J2 ^nQn^n' 

P n,n' 

Qn,n> = <W - KYJl^n'^ir + lp,)U n ^V nfl + 5 <n _ A (r - 7a)^I- A)A KI-a,a]- ( 19 ) 
A 

The indices n, n' are triples of integers, such as (ni, n 2 , n 3 ), labeling the lattice sites, 
while fj, denotes directions, r is the Wilson parameter, K the hopping parameter, U n p, = 
e l 9 aa n p.^ y n * = e «eaA nA ^ w ith q,s the statistical gauge potential and A n ^ the external 
electromagnetic potential. Up is the product of the links surrounding a plaquette, as 
usual, and (3g is related to the three dimensional Us (I) coupling constant by j3a = p^, 
while we denoted by e the four dimensional (dimensionless) electromagnetic coupling of 
C/e(1). In our treatment we will use naive Wilson fermions, so we set r = 0. 

On the other hand, we would like to impose an external homogeneous magnetic field 
in the (missing) x% direction, thus we choose the external gauge potential in such a way 
that the plaquettes in the x±x 2 plane equal B, while all other plaquettes equal zero. This 
may be achieved in several ways, related with each other by gauge transformations. The 
value of B is restricted, due to the fact that the flux through the entire x±x 2 section of the 
lattice should equal an integer multiple of 2tc, so for the x\x 2 plane with area N 2 we have: 
B = j^l, with I some integer between and (The model with integers I between 

and iV 2 is equivalent to the model with integers iV 2 — I which fall into the previous 
category.) We remark here that the maximum value the magnetic filed B can take is n. 
The physical magnetic field B p h ys is related to B through B = ea 2 B p h ys , and the physical 
field may go to infinity, as it should, if the lattice spacing a goes to zero, while B is kept 
constant. 

However, the above considerations are complicated by the fact that the lattice is a 
torus, that is a closed surface, and the Maxwell equation div'B = implies that the 
magnetic flux through the lattice should vanish. This means that, if periodic boundary 
conditions are used for the gauge field,^ the gauge field configuration will have zero total 

8 An alternative way out might be, of course, the use of twisted boundary conditions. 
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flux, so the (positive, say) flux B, penetrating the majority of the plaquettes, will be 
accompanied by a compensating negative flux in a small number of plaquettes. This 
point will become clear in a while. 

To begin with, we consider a configuration giving flux B in each plaquette. In analogy 
with the continuum, one would consider something like: aeA\ = —B112, A 2 = 0, A3 = 0, 
with B = jj^l. However, due to the periodic boundary conditions, the plaquettes in the 
x\x 2 plane starting at the point (ni, n 2 ) with n 2 = N, will not have flux equal to B. To cure 
this, one modifies a little bit the configuration by demanding aeA 2 = B(Nn 1 — y — 
whenever n 2 = N. Thus, all the plaquettes have now flux equal to B, with the exception 
of the one starting at {rii = N, n 2 — N). This plaquette is characterized by flux equal 
to —B(N 2 — 1), which is exactly what is needed to cancel the contribution of all other 
plaquettes. We remind the reader that the value of these plaquettes does not depend on 
£3, so, considering the three-dimensional configuration, we have a constant flux through 
the N 2 — 1 plaquettes and a "Dirac string" of flux of opposite sign through the (N,N) 
plaquette for every X3. 

What we would like to study is the < > condensate, which will depend on the 
lattice point in which it will be calculated. The dependence will be very small if one 
uses sites deep inside the region where the magnetic field is homogeneous. One should 
expect a strong dependence on the site near the location of the string of the large negative 
flux. We have measured the condensate in the middle of the lattice. Another concern is 
the kind of boundary conditions one will use for the fermions. In the time direction we 
used antiperiodic boundary conditions, conforming to common practice. For the spatial 
directions one may use periodic or fixed boundary conditions (the latter meaning that 
the fermion fields vanish on the boundaries). We have used both and found out that the 
results agree for large enough lattices. Finally we chose to work with the fixed boundary 
conditions, so that the fermion matrix feels the least influence from the big negative flux, 
which is located exactly on the boundary. 

4 Lattice Results 

We start with the treatment of the model in the limit (3q — > 00, that is we consider 
the case where the statistical gauge potential is frozen to zero. In this case, the action 
simplifies to: 



S = E ^nQn,n>^n> 




(20) 



where 



aeAi = —Bn 2 , 



aeA 



2 = 5 n2N B(Nn 1 



N N 2 



2"), ^3 = 0, 



2 
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as explained above. We now explain some analytical aspects which may give further 
insight into the model and help us to have an independent check of our results in this 
case. We note that for large enough lattices one may disregard the boundary effects, if 
the condensate is measured away from the regions of inhomogeneity, as explained above. 
Then the gauge field configuration is simplified to: 

aeAi = -Bn 2 , A 2 = 0, A 3 = 0, 

and the action is further reduced to the form: 

S = Y,^nQnn'^n' 

n,n' 

Q nn , = 5 nn , - K 7l (5 n , n+i e"^ - 5 n , t n _ie Wn >) - 

K l2{K',n+2 ~ K',n-2) ~ K lz{K> ,n+3 ~ 5 n>- 3 ) ( 21 ) 

We observe from ( pT|) that in order to take the continuum limit of the lattice action, 
one should consider B on the lattice, such that the maximum BN of Bn 2 should satisfy 
shi-BiV ~ BN, which implies that f\ the magnetic field B should be small enough. If 
we would like a somewhat more quantitative estimate of the maximum allowed magnetic 
field, a trivial calculation of the cyclotron radius can indicate that for B ~ | this radius 
equals two lattice spacings and it decreases for increasing B. Thus the wave functions 
are expected to have structure at distances smaller than one lattice spacing and the 
discretization becomes important. Thus it seems reasonable to stop at this value of B. 

A semi-analytical evaluation of the fermion condensate < ^ (o n2 o)^ (on 2 o) > is achieved 
by first taking the Fourier transform of Q n2 n' 2 for the directions 1 and 3, and then inverting 

n 2 n' 2 ■ 

< *(0n 2 0)*(0n 2 0) >= (75)^2 (22) 

iviv r p ljP3 v 

Qn 2 n' 2 = <5 n2 ^-2iK7isin(pi+Bn 2 )<5„ 2 ^-2iK7 3 sm^^ (23) 

The quantities pi and p 3 in the sum take the values ^rii and ^-(^3 + |) respectively. 
We remark here that this approach parallels the one of Landau when he first derived 
the Landau levels for the motion of a particle in a homogeneous external magnetic field. 
We have checked that the results from the full lattice calculation and the ones from the 
procedure outlined in ( p3|) agree perfectly. 

In figures 4 and 5 we plot A < >= < > |#— < > \b=o versus the magnetic 
field B. In this and the following figures we use the ratio b, defined by b = — . Note 

€.t> \max 

that eB\ max = tc (recall the discussion after equation ([19]) ). 

Figure 4 shows the behaviour of A < >, when the mass parameter m is 0.05 for 
various lattice volumes. We see clearly very strong volume effects, as expected for this 

9 We thank G. Tiktopoulos for a useful discussion on this point. 
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Fermion condensate versus Magnetic field for m=0.05 
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Figure 4: A < > versus the magnetic field for m = 0.05 and various volumes. 

small mass value. As explained in figure 2, the lattice size should be around 100, if one 
is to get strictly zero temperature results. We see that the condensate keeps increasing 
with the volume. The finite temperature effects are stressed emphatically through the 
comparison of the results for the 24 3 and 24 2 x 48 lattices. We should note that we expect 
a non-trivial wave function renormalization of the composite operator < >, which 
can probably account for the difference in the values of the condensate as shown in figures 
1 and 2 on one side and figure 4 on the other. The small B region has qualitatively 
the same behaviour as the corresponding result in figure 1. One may also observe the 
complicated behaviour for large magnetic fields, which is clearly a lattice artifact, as 
explained previously. 

Figure 5 depicts the mass dependence of A < >: Two masses (m = 0.20, 0.10) 
have been used. For m = 0.20 two different volumes have been used (24 2 x 48 and 
36 2 x 72). We recall from the discussion in figure 2 that for lattices of linear extent 25 or 
higher we should get zero temperature results. This is really what is clearly seen in figure 
5: the two lattices, assumed large enough for this mass, give identical results. We have 
also plotted the results for m = 0.10 on the same figure. Two facts have to be pointed 
out: a) In the small B region we see that the condensates for m = 0.10 are bigger than 
the ones for m = 0.20, as anticipated in figure 1. b)The condensates go to zero with the 
mass for fixed volume. This just means that to get meaningful results one should first 
take the thermodynamic limit and after that the massless limit. 
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Figure 5: A < > versus the magnetic field for m — 0.10 and m = 0.20. 



Another issue that has to be settled is gauge invariance, in the sense that the var- 
ious gauge potentials that may be used to generate a homogeneous magnetic field are 
equivalent. In figure 6 we show the result for the condensate for two gauge choices: 
A 1 = -Bn 2 , A 2 = 0, A 3 = 0, and A x = A 2 = A 3 = 0. Lattice volumes 

6 3 and 16 3 have been used. We observe an approximate gauge invariance for the small 
volume and exact gauge invariance for the large volume. The discrepancy in the small 
volume has to do with the fact that only approximate gauge invariance can be achieved 
on small lattices, which are most sensitive to the boundary conditions. 

After the discussion on the free fermion case, we next proceed to incorporate quantum 
fluctuations of the gauge field in the analysis. 

As a first step, we show in figure 7 the condensate for zero magnetic field versus the 
mass m for a 16 3 lattice and three values of the gauge coupling fta '■ 1, 3 and infinite (the 
latter corresponds to the free field). The value 1 is in the strong coupling regime, where a 
condensate will be formed dynamically, while the value 3 is in the weak coupling region, 
where the condensate is small (could be zero if we had dynamical flavors with more than 
four flavours). We observe that the values for (3q = 3 are very close to the ones for the free 
field theory; the extrapolation to zero mass for (3 a = 3 will yield a very small value. On 
the other hand, the condensate for fta — 1 are impressively big, signalling the breakdown 
of the global 577(2) symmetry; the extrapolation to the massless limit certainly points 
to a non-zero value. We should point out here that we should subtract the perturbative 
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Figure 6: < > versus the magnetic field for two gauge choices and two lattice volumes. 



Fermion condensate versus m for zero Magnetic Field, vol=16 A 3 




Figure 7: < > versus m for zero magnetic field and three values of the gauge coupling. 
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Fermion condensate versus Magnetic Field, vol=16 A 3, betag=3. 





Figure 8: < > versus the magnetic field for m = 0.20 and f3 G = 3, for a 16 3 lattice. 



contribution from our lattice result. In that sense our result is preliminary, since it lacks 
the perturbative calculation. We plan to include this in a future publication. 

In the last two figures we depict the results for the full model versus the magnetic 
field. In figure 8 one may see the 16 3 result for (3 G = 3, m = 0.20 These results should be 
compared against the ones in figure 6 for V = 16 3 , Pa — 3, m = 0.20. They are obviously 
very similar, as expected. 

In figure 9 we show the results for 16 3 lattice, fta — 1 and m = 0.20, 0.10, 0.05. We 
observe the same qualitative behaviour with the free case in the small B region. However 
there is a crucial difference: a non-zero condensate is formed, however small the bare 
mass may be, because of the gauge interaction. On the other hand, we find an impressive 
smoothening out of the graph in the large B region; moreover, the difference between the 
results for the various masses decrease with increasing B, so, taken at face value, these 
results imply mass independence of the results in the large B region, as expected on the 
basis of the continuum results. However, in this region we expect many lattice artifacts, as 
explained before, so one should be careful with the interpretation of the subtle continuum 
limit. 
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Fermion condensate versus Magnetic Field, vol=1 6 A 3, Betag=1 . 
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Figure 9: < > versus m for zero magnetic field for m = 0.20, m = 0.10, m = 0.05 
and /9g = 1, for a 16 3 lattice. 
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